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Abstract 

A simple quasi-static model applicable to a wide class of wire media is developed that explains 
strong non-locality in the dielectric response of wire media in clear physical terms of effective 
inductance and capacitance per unit length of a wire. The model is checked against known solutions 
and found to be in excellent agreement with the results obtained by much more sophisticated 
analytical and numerical methods. Special attention is given to suppression of the spatial dispersion 
effects in wire media. 
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I. INTRODUCTION 



Wire media are structured materials formed by many conducting wires embedded in a host 
medium. The wires are normally considered to be very long compared to the wavelength 
in the host medium, but the diameter of the wires is only a small fraction of the lattice 
constant. The known analytical models of wire media 1-11 treat them as crystals of infinitely 
long conducting cylinders. The cylinders may be arranged in different types of lattices 
resulting in different types of anisotropy of the wire crystals. It is known that wire media 
may exhibit strong spatial dispersion, so that the permittivity dyadic e(u, k) in such media 
depends on both frequency and wave vector. For instance, the permittivity dyadic of uniaxial 
wire medium with one set of thin ideally conducting wires oriented along z reads 7 



where ko = uj^e j2 , e and /i are the permittivity and the permeability of the host medium, 
k p is the plasma wavenumber, k z is the wave vector component along z , and I t is the unit 
dyadic in the plane orthogonal to z . 

It is well known that the wire medium supports propagation of transverse electromagnetic 
modes (TEM) which are basically the modes of a multi-wire transmission line. 7,10 Such modes 
propagate along the wires with the velocity equal to the speed of light in the host medium. 
The distribution of the microscopic E and H fields associated with the TEM modes is static- 
like in the planes orthogonal to the wires, with the electric force lines emerging from and 
ending at the surfaces of the wires. It can be easily proven that there is electrical charge 
accumulated on the wires associated with these modes. In Ref. 12 it was shown (for the 
uniaxial wire medium case) that when this charge and the related potential are taken into 
account it is possible to obtain Eq. (1) from simple quasi-static considerations similar to 
those used in Ref. 4. Thus, it was shown that the strong spatial dispersion in wire media 
can be correctly described in a quasi-static approximation. In this paper we extend these 
considerations to a wide class of wire media, and propose an analytical model based on the 
effective inductance and capacitance per unit length of a wire. 

The other motivation for this study is the suppression of the nonlocal effects in wire 
media. In a recent paper by Demetriadou et al. 13 the charge accumulated on the wires 
together with the rather small capacitance of thin wires were identified as the reasons for 
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the spatial dispersion in wire mesh: a metamaterial formed by three sets of wires oriented 
along three Cartesian coordinate axes and joined at the crossing points. A rigorous analytical 
model of such medium was developed in Refs. 9,11. The authors of Ref. 13 make use of this 
model and full wave simulations to justify their main claims. They also propose certain ways 
how to decrease the spatial dispersion effects. The basic idea is to increase the capacitance 
of the wires by periodically loading them with metallic bodies or patches, or alternatively 
to increase the inductance per unit length by coating the wires with a magnetic material. 
Somehow related to this work, it was shown in Refs. 14,15 that for a substrate formed by 
a wire medium slab capped with an array of patches (the so-called mushroom substrate 16 ) 
the response of the wire medium is essentially local. A different strategy to reduce the 
spatial dispersion was reported in Ref. 11, where it was shown that at infrared frequencies 
the plasmonic properties of metals may enable the design of artificial plasmas that mimic 
more closely a continuous local isotropic medium with negative permittivity. 

In this work, we generalize the theories reported in previous studies, 1-11 and propose 
a quasi-static homogenization model that accurately characterizes the nonlocal dielectric 
function of a wide class of wire media (both arrays of parallel wires, and arrays of connected 
wires), including the case where the wires are periodically loaded with conducting metallic 
bodies. In particular, we demonstrate that our analytical theory models accurately the 
electric response of a uniaxial wire medium loaded with patches, and we discuss the physics 
of the suppression of spatial dispersion in such structures. 

II. UNIAXIAL WIRE MEDIUM 

We will start with the simplest possible case of the uniaxial wire medium with one set of 
wires oriented along the z-axis. We will follow the treatment presented in Ref. 12. 

We are interested in the longitudinal (zz) component of the permittivity dyadic. To get 
an expression for it in the quasi-static limit we assume that the radius of the wires r and 
the distance between the wires (the lattice period) a are much less than the wavelength in 
the medium. Let us note that for the model we are going to develop the exact arrangement 
of the wires is not important, it is just enough to know the average distance between a pair 
of neighboring wires in a structure. 

Denoting the average (macroscopic) electric field along z axis in the medium by (E z ), 
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FIG. 1: (Color online) A pair of wires of the uniaxial wire medium without (on the left) and with 
patches (on the right). The integration path used to define Eq. (4) is shown by the blue rectangular 
contour. 

one can write the following relation between this field component and the current in the 
wires l,\ 



{E z ) = {juoL + Z W )I Z + 



dip 

dz ' 



(2) 



where L is the effective inductance per unit length of the wire, Z w is the self-impedance 
of the wire per unit length which accounts for the finite conductivity of metallic wires at 
microwave frequencies or plasmonic behavior at optical frequencies, and <p is the additional 
potential due to charges on the wires. 

This relation can be obtained integrating the microscopic electric field over a path shown 
in Fig. 1. The path goes first along the surface of a wire then to the middle line in a pair 
of two neighboring wires, then along this middle line and, finally, back to the surface of the 
wire. The circulation of the microscopic electric field E(x, z) over this path reads 

z+Az z+Az a/2 a/2 

(j)E d\= J E z (r ,z')dz' - J E z (a/2, z') dz' + J E x (x, z + Az) dx - J E x (x,z)dx. 



(3) 

The first integral in this relation represents the voltage drop along the surface of the wire 
and, therefore, can be expressed in terms of the wire current and the wire self-impedance 
per unit lenght. The second integral is the voltage drop along the symmetry line shown in 



Fig. 1. In the same manner as it was done in Ref. 4 we relate this voltage drop with the 
macroscopic electric field in the medium. After doing this the circulation of the electric field 
reads (when Az is small enough) 

a/2 

<j> E • dl = (Z W I Z - (E z ))Az + <p(z + Az) - tp(z), where <p(z) = J E x {x, z) dx. (4) 



ro 



The electric field circulation equals minus the time derivative of the magnetic flux that 
penetrates the area bounded by the integration path: § E • dl = — juoQ = —jujLI z Az, from 
which we immediately get (2) when Az — > 0. 

In general, the effective inductance L depends on the specific microstructure of the system 
(e.g. if the wires are coated or not with some material). In the particular case in which 
the wires are conducting cylinders (with no material coating), it was shown, 4 by calculating 
the magnetic flux of a pair of neighboring wires in the quasi-static approximation, that L 
verifies: 

L = ^log / - . (5) 

2tt S 4r (a-r ) V ; 

It may be verified that the above formula also applies to the case where the wires are loaded 
with metallic patches (Fig. 1, right). 

The additional potential caused by the charges on the wires can be found by placing 
a linear charge density p on the wires and by calculating the corresponding electrostatic 
potential ip created by the fluctuating part of the microscopic electric field. Thus p is 
responsible for the electric field component orthogonal to the wires. We introduce an effective 
capacitance C per unit length, such that it verifies: 

= —■ (6) 

Notice that the considered capacitance is calculated by placing an identical linear charge 
density over the wires (differently from the traditional definition of capacitance, which as- 
sumes that charge density over two conductors is antisymmetric). In the same manner as 
the inductance, the capacitance depends on the microstructure of the system. In the quasi- 
static limit a pair of charged wires (with no attached conducting bodies) induces the field 
(see Fig. 1) 

(7) 



E x = — — — 

2ne Vl 



x a — x 
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This expression has the same form as the one used in Ref. 4 for the quasi-static magnetic 
field of a pair of lines of current. Therefore, for this particular case the capacitance is given 
by 

11, « 2 

°S ~r~~r ~~t- (8) 



C 2ne & 4r (a-r )' 
The capacitance for a system of wires loaded with conducting patches (Fig. 1, right) is 
calculated in Appendix A. 

Considering now a monochromatic plane wave of current excited in the crystal, the cur- 
rents in the wires can be written in the form 

I g (z)=I e-* k ' z , (9) 

and thus the linear density of the charge associated with the currents verifies 

P (z) = -L d JM = h Iz{z y (10 ) 
ju) dz uj 

These charges are responsible for the electric field component orthogonal to the wires. 

Hence, the relation (2) can be rewritten in terms of the effective inductance and of the 
effective capacitance per unit length of the wire as 

(E z )= (jujL + Z w + -£-)l g . (11) 

Already in this expression one can identify the spatial dispersion term proportional to the 
square of the z-component of the wave vector. 

The macroscopic polarization current in wire media is the average of the currents in 
separate wires. Let A cc \\ be the average area in the xy plane per one wire of the crystal. 
Then the macroscopic polarization current is J z = I z /A cc \\. The macroscopic displacement 
field is D z = Eq (E z ) + J z /{juj). Therefore, after some algebra we find that the longitudinal 
component of the permittivity dyadic is given by 

k 2 

ffi = 1 _E do) 

eo *8-#A*)-*2/n 2 ' 1 ) 

where k"i = (j,o/(A cc \\L), n 2 = LC/(eo/J>o), £ = (Z w /L)y/eojIo. It may be easily checked 
that the above formula reduces to Eq. (1) in the case of perfectly conducting straight wires 
(Z w = 0) [also, for unloaded wires n = 1 as is seen from Eqs. (5) and (8)]. More generally, 
when the wires are characterized by the complex permittivity eoe m (e.g., thin plasmonic 
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rods at optical frequencies), the impedance Z w is given by, 

Z w = - 2 7 77 > (13) 

where r is the radius of the rods. It may be easily verified that in this scenario Eq. (12) 
reduces to formula (16) of Ref. 10, which was calculated using a local field based approach. 
Thus, Eq. (12) generalizes the previous homogenization models of the uniaxial wire medium. 

Nevertheless, it is worth noting that the expression for the plasma wavenumber obtained 
in the present paper differs from the one derived in previous works. 3 ' 10 Namely, under the 
approach developed above we have 

(W 2 = : • (14) 

|qct - 

° 4ro(a— ro) 

In Refs. 3,10, under a thin wire approximation, it was obtained that 

2-7T 

(M 2 « 7T- (15) 

1 p ; 0.5275 + log ^ 1 } 

One can notice that (15) gives unphysical results for any r /a > (27r)~ 1 exp(0.5275) ~ 
0.27. Contrary, Eq. (14) gives a physically sound result in the limit ro — > a/2 when the 
surfaces of two wires touch: It predicts an infinite growth in the magnitude of k p in this 
limit. It can be also checked numerically that the accuracy of (14) is better than (15) when 
r PS 0.1a or larger, whereas the opposite behavior is observed for r < 0.05a. Nevertheless, 
both formulas have the same asymptotic behavior when r — > 0. At r /a = 0.05 (this ratio 
has been used in our numerical simulations that are discussed in Section IV) the formulas 
(14) and (15) overestimate the plasma frequency by about 3%. 

Another asymptotic expression for the normalized plasma frequency which is often cited 
was obtained in Refs. 1,2 but even for rather small wire radii its accuracy is worse than that 
of (14) and (15). Also, it does not predict the infinite growth of k p when r — > a/2. 

It should be emphasized that Eq. (12) is in principle valid for a wide class of wire media 
(e.g. wires with attached conducting bodies). The parameters C and L depend on the 
specific microstructure of the system. The magnitude of the spatial dispersion term k^/n 2 



in (12) can be reduced by increasing the value of n — a/ LC /(eqUo). This quantity has 
the meaning of slow-wave factor for quasi-TEM waves propagating along the wires. As 
mentioned before, for unloaded straight wires n — 1. As discussed in Ref. 13, the capacitance 
C can be increased by loading wires with metallic patches and the inductance L can be 



increased by placing wires in ferromagnetic shields. An alternative way to increase the 
inductance is to use helices instead of straight wires. Associated bi-anisotropy in helix 
medium can be compensated if both right- and left-handed helices are used. 

Attaching metallic or dielectric bodies to the wires also changes the transversal compo- 
nents of the permittivity dyadic. We will study this effect with more details in Section IV. 



III. WIRE MESH 

The (3D) wire mesh is a wire crystal formed by three mutually orthogonal sets of wires 
joined at the intersection points. The electromagnetics of such metamaterial have been 
studied in several recent works. 9 ' 11 ' 17,18 In the following derivation we assume a cubic lattice, 
but after a straightforward generalization the same method can be applied to structures of 
more complex geometries. Similar to the case studied in section II, metallic or dielectric 
bodies may be attached to the wires. 

In the wire mesh we get three components of the polarization current related with the 
currents in three orthogonal sets of wires. The currents in the wires are related to the 
average electric field in the medium in a manner similar to the uniaxial case: 

(Ex) = {jujL + Z W )I X + ^, (16) 

(Ey) = (juL + Z W )ly + ^ : (17) 

(E z ) = ( JU jL + Z w )I z + ^. (18) 

Because the wires are joined at the crossing points they are locally under the same potential, 

that is why we have the same in all three equations. But the currents in three sets of 

wires can differ and that is taken into account by the variables I x , I y , and l z . 

Let us consider a unit cell of the wire mesh with three intersecting connected wires. The 

total charge q accumulated on these three wires per unit cell can be found as 

a f dl x ^ dl y ^ dl z \ 
ju \ dx dy dz J 

Because the wires are electrically connected and their effective capacitance per unit length is 

the same, this charge is equally distributed among the three wires in the unit cell. Therefore, 

for the linear charge densities on the wires we have in a vicinity of the unit cell 

q If dl x dl v dl z \ 

Px = Py = Pz = Ta = -^\-d^ + lty- + lu)- (20) 



Using the same notation for the effective capacitance of a wire as above we can write the 
potential ip as 



1 



<P = 



dl x dl y dl z 



1 



3juC \ dx dy dz J 3uC 
where we have taken into account that the currents on the wires change on average as 

I n = I°e-i k " n , n = x,y,z. 



(21) 



(22) 



Now we can substitute this expression for the additional potential into (16)-(18). Doing 
this we obtain the following system of equations: 



(E x ) = (j U ;L + Z w + 7 ^)I x + 7 ^(k y I y + k z I z ), 



(E y ) = (juL + Z w + 



3jujC 3juC 

k, 



ly ~t~ q . n {kzlz + k x I x ), 



3jcuC J y 3juC 
(E z ) = (jcoL + Z w + -^-)I z + -^-(k x I x + k y I y ). 



(23) 
(24) 
(25) 



3juC 3juC 

By introducing a vector of currents I = J x x + I y yo + I z z we rewrite this system in a more 
compact form using dyadics: 

kk 



(E) = 



(juL + Z W )I 



3juC 



I, 



(26) 



where I is the unit dyadic and kk = k <8> k is the dyadic (tensor) product of two vectors. 
Now it is only a matter of inverting the dyadic in brackets of (26) to get the permittivity 
dyadic of the wire mesh. 

The average polarization in the medium is P = 1/ (jujA cc \\) + P t , where P t accounts for 
additional polarization due to finite thickness of the wires or metallic bodies attached to 
the wires. For a crystal of cubic symmetry we can write P t = s (e t — 1) (E), therefore the 
displacement vector D = EoE t (E) +I/(juA cell ), and 



e(uj,k) 



(juL + Z w )I + 



kk 

3juC 



-i 



(27) 



or, after some dyadic algebra, 
f (w, k) 



k 2 



k 2 -jtk J 3n 2 [e -^k ][k 2 -^k - A;V(3n 2 )]' 



(28) 
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where we use the same notations as in (12), and k 2 = k 2 + k 2 + k 2 z . The obtained permittivity 
dyadic can be also rewritten as 

e(c<j,k) . . /= kk\ . , x kk . . 

} =e tr (tu)(l- — )+eUcu,k)—, (29) 



where 

£tr(w) = 5t " (30) 

£lo( ^ fc) = £t - ^-^-A:V(3n 2 ) - (31) 
It can be verified that for the mesh of thin plasmonic rods without loading [for which Z w 
is given by Eq. (13)], the relations (30)-(31) transform to the ones presented in Ref. 11 with 
the parameters e t = 1, k p = (3 P , and identifying the numerical coefficient lo from the same 
reference with L = 3n 2 . 



IV. UNIAXIAL WIRE MEDIUM LOADED WITH PATCHES AND SUPPRES- 
SION OF SPATIAL DISPERSION 

Recently 13 it was proposed to load the wire mesh with metal patches to increase the 
effective capacitance of the wires per unit length and decrease the related spatial dispersion 
effects. This proposal was supported by numerical simulations. Here, we will apply our 
general analytical model to the particular case of a uniaxial wire medium loaded with metal 
patches. For this purpose we just need to determine what is the effective capacitance C 
introduced in Section II in the presence of patches. The details of calculation of this capaci- 
tance are described in Appendix A. Here we give the result: C = C wirc + C pat ch, where C w i rc 
is the wire capacitance given by (8) and 

Cpatch = hhg(Z%y (32) 

where w is the width of the square patches periodically attached to the wires and separated 
by the distance h along z, and d = a — w is the gap between two adjacent patches on a 
pair of neighboring wires. Thus, the permittivity dyadic of the uniaxial wire medium loaded 
with patches is given by 

10 



where we keep the same notations as in Section II. The transverse permittivity e t is mostly 
determined by the patches when w ^> r and it can be found as the permittivity of a stack 
of capacitive grids separated by h one from another. With the help of the known theory of 
such grids 19 it can be found that 

El = l + _log^c-j. (34) 

The accuracy of (32) and (34) is better for small gaps and for large values of h/a. 

In the limit d — > the effective capacitance behaves as C ~ l% ^f and, therefore, can 
be arbitrarily large if the gap between two adjacent patches is made small enough. On the 
other hand, the transverse permittivity e t grows under the same limit as e t ~ ^ log (%)■ 
The square of the slow-wave factor n 2 is proportional to the effective capacitance, therefore, 
by increasing the width of the patches one can discard the spatial dispersion term in the 
right-hand side of (33) while keeping e t at a reasonable level (this is possible because e t 
grows more slowly when d — > 0). An explicit expression for the slow- wave factor under the 
mentioned limit reads 

9 LC „ 16w 

n = — ~ 1 + m ^2 - 35 

e /2 nh{k v dy 

In fact, we have numerically checked that this simple expression works quite well for gaps 
of width d < 0.2a. 

To illustrate the suppression of the spatial dispersion in the considered wire media, we 
have calculated the dispersion diagrams for several configurations using our quasi-static 
model, the transfer matrix method described in Appendix B, and the eigenmode solver of 
CST Microwave Studio. The structure was assumed lossless in the simulations (all metallic 
components are perfectly conducting so that Z w = 0). The transfer matrix formalism 
developed in Appendix B is based on the assumption that in between two patch grids the 
electric field is a superposition of TEM and TM modes. 7 The fields on the interfaces of 
each patch grid are linked by a grid impedance and by an additional boundary condition, 20 
consistent with the formalism described in Refs. 14,15. The obtained results are presented 
in Fig. 2. 

In Fig. 2(a) and Fig. 2(c) the dispersion diagrams obtained from the quasi-static model 
and the numerical simulations are shown for a set of the propagation angles with respect to 
the axis of the structure: a — 0, 30°, 60° [for the other parameters of the structure refer to 
Fig. 1; in these plots the wave vector is k = k (sino;x + cosaz )]. The dispersion curves 

11 




FIG. 2: (Color online) Dispersion diagrams for a uniaxial wire medium loaded with patches ob- 
tained using two analytical models and numerical simulations for different propagation angles a 
with respect to the z-axis. Only the branches associated with the quasi-TEM and TM modes 
are shown. Panels (a) and (c): quasi-static model vs. numerical simulations: (a) w = 0.5a, (c) 
w = 0.9a. Panels (b) and (d): transfer matrix model vs. numerical simulations: (b) w = 0.5a, 
(d) w = 0.9a. On all 4 panels the solid lines represent the analytical results and the symbols 
correspond to the results of numerical simulations; the values of the propagation angles are coded 
in color: a = 0: blue lines and circles; a = 30°: magenta lines and triangles; a = 60°: red lines 
and crosses. The other parameters in all 4 cases: ro = 0.05a, h = a. 

predicted by the quasi-static model are depicted with solid lines while the results of the 
numerical simulations are represented by symbols. In the example of Fig. 2(a) the patch 
width has been set equal to w = 0.5a, while in Fig. 2(c) the patch width is w = 0.9a. In 
both cases the theory and the simulations predict the existence of two dispersion branches 
associated with extraordinary waves, i.e., with the quasi-TEM and TM modes, as well as a 
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dispersion branch associated with the ordinary (TE) wave whose dispersion is not depicted 
in Fig. 2 (there are also other higher order modes at higher frequencies, but we are not 
interested in them). We call the high-frequency branch "the plasmon mode" because for 
a = this branch corresponds to the longitudinal plasmon-type wave propagating along 
the axis of the structure. On the other hand, the low-frequency branch for a = belongs 
to an ordinary transverse wave which is not affected by the wires (but it is affected by the 
transverse permittivity e t of the medium). 

From Fig. 2(a) one can see that for the moderate-size patches the quasi-static model 
works surprisingly well even when ka approaches n. The small difference in the frequencies 
of the plasmon-type modes predicted by the theory and the simulations at ka = is due to 
the asymptotic nature of the formula for the plasma wavenumber that we use (the discussion 
on this is given in Section II). For larger patches (Fig. 2(c)) the quasi-static model does not 
predict appearance of a band gap at a = and ka = ir. This is expected since in the model 
the capacitive loading on the wires is assumed to be effectively uniform along the wires. 

Fig. 2(b) and Fig. 2(d) display the same dispersion diagrams but with the quasi-static 
model replaced by the transfer matrix model described in Appendix B. One can see that 
this model wrongly predicts a completely flat dispersion for the plasmon mode propagating 
along the z axis (a = 0), independently of the patch size. This is in disagreement with 
the numerical simulations, as is seen from Fig. 2(b). Indeed, the formalism developed in 
Refs. 14,15 is only valid when the gap between the patches is small, because otherwise other 
higher modes can be excited near the connections of the wires to the patch grid, and in 
such conditions it is not possible to consider that the microscopic field in the vicinity of the 
connection points are a superposition of TM and TEM modes of the unloaded wire medium, 
as assumed in Refs. 14,15. Consistent with this observation, it is seen in Fig. 2(d), that for 
larger patches and (or) larger angles of propagation the disagreement is less pronounced. 
Another characteristic feature of the transfer matrix model is that it is able to predict the 
existence of the above-mentioned bandgap. This is because the transfer matrix model takes 
into account the granularity of the structure along the z axis. 

The suppression of the spatial dispersion effects is evident if we compare Fig. 2(a) with 
Fig. 2(c). Indeed, the latter case corresponds to a larger patch width [w = 0.9a), and 
consequently the slope of the dispersion curve associated with the longitudinal mode (the 
plasmon mode at a — 0) is very small. To justify this effect and also to check the accuracy 
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FIG. 3: (Color online) The square of the slow- wave factor as a function of a/d (logarithmic scale). 
The lines represent the result of the quasi-static model, the symbols correspond to the values of 
n 2 extracted from the numerical simulations. Blue dotted line and crosses: ro = 0.05a, h = a/3; 
red solid line and circles: ro = 0.05a, h = a. 

of the quasi-static model near the origin of the Brillouin zone for a wide range of values 
of the gap, we have extracted the values of the slow wave factor n from the results of the 
numerical simulations slightly above the point ka = and compared them with the value of 
n given by the analytical model. The results of this extraction are presented in Fig. 3. From 
this figure we see that despite its simplicity, the quasi-static model predicts very well the 
trend in the growth of n 2 when the the gap between the patches decreases. The agreement 
tends to improve for larger values of h/a. 

V. CONCLUSIONS 

In this paper we have developed a quasi-static analytical model of wire media applicable 
to a wide class of structures, and in particular we have considered uniaxial and isotropic 
wire crystals, which may be loaded with metallic patches. Because the developed model is 
defined in simple physical terms of the effective inductance and capacitance per unit length 
of a wire it can be readily extended to other wire structures of more complex geometries. 
The model accounts for the finite conductivity of the wires so that it can be applied when 
the metallic wires become plasmonic (consistent with the results reported in Refs. 10,11) or 
when the wires are uniformly loaded with arbitrary complex impedances. In particular, we 
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have studied with details the electrodynamics of uniaxial wire media loaded with patches, 
and demonstrated with full wave simulations that the proposed quasi-static model describes 
accurately the properties of the system in the long wavelength limit. Consistent with the 
analysis of Refs. 13-15, it was shown that the presence of the patches may result in a dramatic 
reduction of the nonlocal effects. For the case of unloaded wire media, we have demonstrated 
that the quasi-static model yields the same expressions for the dielectric permittivity tensors 
as those obtained by much more sophisticated methods. 5 ' 6,8 ~ n Thus, we have proven that the 
strong spatial dispersion in wire media is a quasi-static effect. Although this fact has already 
been noticed, 12 the presented research extends the results obtained in Ref. 12 and allows 
for analytical and quantitative studies of the possibilities to control the spatial dispersion in 
wire media. 
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APPENDIX A: 



As is seen from Fig. 1 (right) that depicts the path along which we calculate the circulation 
of the electric field, the capacitance in question can be calculated if we find the electric field 
in the region close to the gap between two patches on the neighboring wires. Indeed, the 
circulation integral (3) in the presence of patches has to be modified as follows 

z +h z +h a/2 a/2 

<j) E ■ dl = J E z (r , z) dz — J E z (a/2, z) dz + J E x (x, z + h) dx — J E x (x,z )dx, 

zq zo w/2 w/2 

(Al) 

where zq is at the location of an arbitrary plane of patches. We choose the integration path 
so that it first goes along the surfaces of the wire and the patch till the gap, then across the 
gap till the symmetry line (Fig. 1), then along that line till the second gap and then across 
this gap back to the patch and the wire. 

One can see that the first integral in the right-hand side of (Al) still represents the same 
quantity as in the unloaded uniaxial wire medium and is related to the finite conductivity 
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of the wire. The integrals along the surfaces of the patches are not shown in (Al) as 
in the following we consider the patches to be ideally conducting. This may be a good 
approximation because in practice the wire impedance dominates. We can express the 
second integral of (Al) as 

zo+h 

J E z (a/2, z) dz = (E z )h - ip z {z + h) + ip z (z ), (A2) 

ZQ 

where the two last terms account for the strong non-uniformity of the z-component of the 
microscopic electric field in the vicinity of the gaps and can be defined (in the unit cell 

z < z < z + h) as 

zu+h/2 

tp z {z) = J (E z (a/2, z) - (E z )) dz. (A3) 

z 

Notice that when there are no patches the microscopic field changes smoothly along z; that 
is why in Section II we could simply relate the second integral of (Al) with the macroscopic 
electric field. The two remaining integrals can be written in terms of 

a/2 

<Px(z) = J E x (x,z)dx. (A4) 

w/2 

Substituting the above expressions for the integrals into (Al) and comparing it with (4) 
with Az = h we see that the additional potential at the plane z = z has to be 

<p(zo) = <Px(zo) + <Pz(zo)- (A5) 

At this point it is worth reminding that the additional potential as we define it and 
use it in Sections II and III is essentially a macroscopic quantity: It changes slowly and 
smoothly along the wires. Therefore, Eq. (A5) can be understood as the definition of the 
averaging procedure for the additional electric field (represented by both non-uniform E x 
and E z components) that appears because of the periodical non-uniformity in the charge 
distribution introduced by the patches. 

It is clear that for wide patches this additional potential is mainly determined by the 
fluctuating part of the microscopic field in the vicinity of the gap at x = a/2, z = z . Thus, 
to simplify the problem we may first neglect the effect of the charges sitting on the wires 
(nevertheless, we will later add a correction term taking the wires into account). Second, 
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because tp(zo) depends only on the non-uniform part of the field we may neglect the effect 
of all other planes of charged patches except the plane z = zq. We may do so because when 
h ^> d = a — w the field produced by the other planes of patches is practically uniform in 
the vicinity of the gap we are interested in. Therefore in the following we consider only a 
single array of charged patches and discard all the patches that are not co-planar with the 
patch at z = z . 

In what follows, we will solve the enunciated electrostatic problem and calculate the 
effective capacitance per unit length of a wire with patches. In order to obtain a closed-form 
analytical solution we will make an additional simplification: we replace all the patches 
centered at the same ^-coordinate with a single metal strip of the same width. Thus, 
we obtain a grid of metallic strips (geometry of the problem becomes invariant along y) 
separated by the same gap as the array of patches. 
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FIG. 4: (Color online) The original domain (a) and the domain obtained after the conformal 
mapping (b) defined by Eq. (A6). The solid red lines represent the perfect electric conductor 
boundary (PEC) (surface of the patch), the dashed blue lines represent the perfect magnetic 
conductor (PMC) boundaries that impose necessary symmetries. The point vq corresponds to the 
origin of the domain (a). 



Taking into account the symmetry of the excitation and the periodicity of the grid we 
arrive at the two-dimensional problem shown in Fig. 4(a). In this figure we define the local 
coordinate system follows. The middle point of the gap is at x' = z' = 0. The 

patch is modeled as an infinitely thin PEC (perfect electric conductor) strip, which starts at 
x' — d/2 and continues to the point x' = a/2 which is at the middle line of the patch. The 
strip is charged. The PMC (perfect magnetic conductor) boundaries shown in the figure 
enforce the symmetries mentioned above. 
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The distribution of the electric potential in this system can be found with the conformal 
mapping approach. 21 One can verify that the following analytical function of the complex 
variable Z = x' + jz' 



C(Z) = u(x'i z') + jv(x', z') = cos 



-i 



2 cos ^ - cos 22 + 1 

a a 

COS ^ + 1 



(A6) 



maps the domain shown in Fig. 4(a) into the domain of Fig. 4(b) in which the solution for 
the electric field is trivial. One can see that the curves v(x f , z') = const are the equipotential 
contours and the curves u(x', z') = const are the force lines of the electric field of the 
problem. It can be checked that at the the patch the potential defined in this way vanishes: 
v(x',0) = for d/2 < x' < a/2. 

Therefore, the voltage drop between a distant point on the z' axis (which is a point on 
the integration path shown in Fig. 1) and the edge of the patch is 



v (0, z') = cosh 1 



2 cosh - cos — + 1 



1 + cos 



ird 



(A7) 



When z' >> d this voltage asymptotically behaves as 



,, , / nd\ 2nz' 

v(0, z) ~ 2 log sec — M . 

\ 2a J a 



(A8) 



The linearly growing term of (A8) corresponds to a uniform electric field far away from the 
gap. Such a smooth field is already taken into account by the first term of (A2). Therefore, 
the additional potential we are looking for must be defined as 



if(z ) = lim ( v(0,z') 



2nz' 



2 log I sec — | . 

2a J 



(A9) 



On the other hand, the total charge per unit length of the strip is given by 



Q = Ae [u(a/2, 0) - u(d/2, 0)] = Ans , 



(A10) 



where the coefficient 4 accounts for the fact that the domain of Fig. 4(a) includes only 1/4 
of the total surface of the strip. Hence, the effective capacitance per unit length of a strip is 
Cstrip = Q/f{zo) and the effective patch capacitance per unit length of a wire with patches 
is Cp atc h = (w/h)C stT ip, which is given by Eq. (32). The total capacitance per unit length 
of a wire with patches is approximated as a sum of the wire capacitance C wire from (8) and 

Cpatch • 
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APPENDIX B: 



The uniaxial wire medium periodically loaded with patches may be regarded as a layered 
structure. Thus, it is possible to apply the standard transfer matrix method to it, provided 
we are able to characterize the fields in one cell. As demonstrated next, this can be done by 
generalizing the formalism developed in Refs. 14,15, which is based on the assumption that 
in the regions in between two arbitrary adjacent patch grids the microscopic fields can be 
written as a superposition of TM and TEM modes of the unloaded wire medium. 

Suppose that the wires are directed along z and that the magnetic field is along the y- 
direction (the TM polarization). The electric field components are E x and E z and the wave 
vector is k = k x x + k z z . The fields in the region 0<z<h(z = 0at the patch array) in 
between two arrays of metallic patches can be decomposed into four waves: 

n„ H - A + P "fTM2 I 4~ p +7tm2 _i_ D+ p "7temz I D- „+7temz mi \ 

l IO rl y — ^TM 6 ' ^TM 6 ' TEM e ' TEM > K 01 ) 

E x = ^L-± (voHy) = -^L- [ 7TM (^TM e_7 ™ Z - ^tm^ 7 ™*) + 

7tem (^ EM e- 7TEMZ + ^ EM e +7TEM2 )] , (B2) 

—k 

^ z = C TM U (^TM e 7 ™ 2 + ^TM e+7 ™ Z ) > (B3) 
S zz K Q 



where 7tm = \Jkl + kl~ fcg, 7 TEM = jk , and e™ = e k 2 x /(kl + k 2 x ). It may be easily 
recognized that the above expressions correspond to a superposition of the standard TEM 
and TM modes supported by the unloaded wire medium. 7 

The tangential components of the fields at the planes z = h~ and z = h + are linked by 
the transfer matrix of the patch array: 



E x 




(B4) 

=h+ \ "° / \ ° / z=h- 



where y g = 2je ^ log (esc (ff)) is the normalized effective admittance of the patch 
array 14,15 On the other hand, it was shown in Ref. 20 that since the microscopic surface 
charge density must vanish at the connections between the wires and the grid, the following 
additional boundary condition (ABC) must be verified: 

ko£o ~Jz~ + kxVo ~Iz~ = °' Z = ° + ' Z = h ~ : 
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where i] is the free-space impedance. Using this ABC in Eqs. (B1)-(B3), it is possible to 
obtain the coefficients associated with the TEM wave (B^ EM ) as a function of the coefficients 
associated with the TM wave (A^ M ). Then, A^ M can be expressed in terms of the tangential 
electric and magnetic fields at the z = + plane. Proceeding in this manner, it is possible to 
obtain after lengthy but straightforward calculations the following transfer matrix relation 
for the layer of the wire medium in between two patch grids: 

E x \ I E, r 



M 



z=h 



\VoHy 



(B6) 



z=0+ 



M 



where the matrix M is 

m 2 i mn 

with the elements given by the following formulas: 

(e - e™) 7tm sinh (7tm^) cosh (7tem^) + £™7tem cosh (ttm^) sinh (7tem^) 



(B7) 



mn 



m i2 = 



rM sinh (7tm^) cosh (7tem^) + gjj , 
(e - e™) 7tm sinh (7™^) + £™7tem sinh (7temM 

_1 J7tem7tm sinh (7tm^) sinh 1 

k (e„ - e™) 7tm sinh (7tm^) + £™7i 



m 2 i = -jk Q 



2{e 



_TM\ r TM 



° zz j ° zz 



tem7tm sinh (7tm^) sinh (7tem^) 

7tm sinh (7tm^) + £™7tem sinh (7tem^) ' 

-1 + cosh (7™^) cosh (7tem^)] 



(B8) 
(B9) 



(e - e™) 7tm sinh (7tm^) + £™1tem sinh (ttem^) 

sinh ( lTEM h) sinh ( lTM h) Ueo - e™) 2 + (e™) 



7TEM 
7TM 



(e - e™) 7tm sinh (7tm^) + £™7tem sinh (7temM 



(BIO) 



It may be verified that det(M) = 1. From Eq. (B6) and Eq. (B4) it is clear that the global 
transfer matrix of the system is: 



M 



• M. 



(Bll) 



1 

~Vg 1 

As is well known, the dispersion characteristic of the Bloch-Floquet modes supported by the 
periodic structure verifies cos (k z h) = tr(M g ) /2, where tr(...) represents the trace of the 
matrix. Thus, it follows that dispersion equation for the Bloch-Floquet modes is: 



v 

cos (k z h) = mn - 7T m i2> 



(B12) 
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The transfer matrix dispersion diagrams calculated section IV were obtained using the above 
equation. 
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